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We consider the possibility that the DAMA signal arises from channeled events in simple models 
where the dark matter interaction with nuclei is suppressed at small momenta. As with the standard 
WIMP, these models have two parameters (the dark matter mass and the size of the cross-section), 
without the need to introduce an additional energy threshold type of parameter. We find that they 
can be consistent with channeling fractions as low as about ~ 15%, so long as at least ~70% of the 
nuclear recoil energy for channeled events is deposited electronically. Given that there are reasons 
not to expect very large channeling fractions, these scenarios make the channeling explanation of 
DAMA much more compelling. 



PACS numbers: 95.35. +d 

I. INTRODUCTION 

The annual modulation signal seen by the Dark Matter 
(DAMA) collaboration, now significant at the 8a level 
[TJE], provides an intriguing puzzle. If one assumes that 
this signature is due to a standard weakly interacting 
massive particle (WIMP) elastically scattering off nuclei, 
then it follows that other direct detection experiments 
should also have seen large numbers of signal events. This 
is, however, not the case. As a result, one is left to explore 
dark matter scenarios going beyond the standard elastic 
WIMP paradigm. Indeed, it is not currently possible 
to rule out a dark matter explanation for DAMA in a 
model independent way. The list of known possibilities 
still consistent with all available data, however, is in fact 
rather limited. 

One scenario which has received a great deal of atten- 
tion is that of inelastic dark matter (3J. In this picture, 
one assumes that in each scattering event the dark matter 
particle makes a transition to a higher energy state, sep- 
arated from the original by a splitting S. This splitting 
effectively removes scattering events at sufficiently small 
nuclear recoil energies, where standard WIMP event rates 
become large and cause the greatest tension with data. 
Another possibility is that the events at small recoil en- 
ergies become suppressed by a dynamical form factor 
[U [5] . Yet another option is that the scattering proceeds 
through a resonance [6]. 

In all of the above cases, however, great care must be 
taken to ensure that the DAMA signal is present, with- 
out resulting in an over-abundance of events at other 
experiments. In inelastic dark matter, a splitting is 
added solely to remove events at experiments with light 
nuclei. This splitting cannot be too small, lest the low 
energy suppression be insufficient, but it also cannot be 
too large, or else the DAMA signal disappears. Hence, 
for any given dark matter mass, 6 must fall in a nar- 
row range, and additional model building is required to 
motivate the coincidence associated with its size. 

In the dynamical form factor scenario, careful model 
building must be undertaken in order to ensure that the 



required low energy suppression is sufficiently rapid. In 
explicit models constructed so far, here too there is a 
parameter (analogous to 5) that must be carefully cho- 
sen to make the DAMA signal consistent with the null 
experiments. 

While such model building avenues are currently still 
viable for explaining DAMA, perhaps a nicer resolution 
would be obtained if there were a reason found for DAMA 
simply being the most sensitive direct detection experi- 
ment at low recoil energies, where WIMP signals tend 
to be peaked. In fact, "channeling" has already been 
proposed as an explanation of this type |5J M DZ3 ITT] . 
Most direct detection experiments, including DAMA, do 
not measure the complete recoil energy of a nucleus in a 
scattering event, but only some fraction of it. In particu- 
lar, DAMA measures only the fraction deposited electro- 
magnetically. The complete recoil energy is then inferred 
through multiplication by an experimentally determined 
"quenching factor" . DAMA, modeling their quenching 
factor to be the same for all recoiling Iodine nuclei, de- 
termines it to have a value of 0.09. The idea behind 
channeling is that, due to the crystal structure of the de- 
tector, some fraction of low energy events may actually 
have a quenching factor closer to one. This is the expec- 
tation for nuclear recoils which travel sufficiently along 
an axis of the crystal. For such "channeled" events, a nu- 
clear recoil which was thought to have had an energy of, 
say <~ 20keV, actually had an energy closer to ~ 2keV. 

In this letter, we will reconsider the possibility of ex- 
plaining the DAMA data through the assumption that 
some fraction of recoiling nuclei are being channeled. 
Since we believe that current theoretical estimates for 
the amount of channeling are most likely not reliable, our 
goal will be to try to explain the DAMA data using as 
little channeling as possible. We will show that unlike the 
case of a standard WIMP, the simplest momentum de- 
pendent couplings can allow for an appreciable decrease 
in the amount of channeling required to explain the data, 
to as little as ~ 10%. In addition, for channeling frac- 
tions of at least ~ 15%, these couplings can allow for 
channeled events to have a quenching factor as small as 
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~ 70% and still satisfy all constraints. 

No complicated model building is needed to arrange 
for the required momentum dependence. For example, 
it is sufficient to have a neutral dark matter particle, X, 
whose leading interaction with a new GeV mass gauge 
boson A M therefore becomes the lowest dimension gauge- 
invariant interaction C — j S F flv d IJ 'X*d v X. Kinetic mix- 
ing between and hypercharge then couples X to nu- 
clei in a momentum dependent fashion. The mass of the 
gauge boson and the scale A ~GeV are both absorbed 
into the overall scattering cross-section, and so consti- 
tute a single parameter. 

The scenario discussed here therefore offers a solution 
to the DAMA puzzle that does not require extra energy 
threshold parameters, and is viable with more conserva- 
tive and realistic versions of channeling. Finally, in the 
event that heavy element experiments such as CRESST 
[H EES], XENON PJ], and KIMS |Hj fail to find evi- 
dence for dark matter in the near future, this scenario 
may remain as an attractive viable option. 



II. REVIEW 

The event rate per unit detector mass per unit recoil 
energy at a dark matter direct detection experiment is 
given by (see e.g. [HUE]) 
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In this expression, Nt is the number of target nuclei per 
unit mass, p x is the local dark matter density and tudm 
is the mass of the dark matter particle. The integration 
is over the distribution of dark matter velocities f(~v) 
relative to the Earth; the lower bound v m i n is given by 
the minimum velocity a dark matter particle must have 
in order to cause a nuclear recoil with energy Er. In 

particular, we have v min = ~\ J&*^ , where tojv is the 
mass of the recoiling nucleus, and /i is the dark matter- 
nucleus reduced mass. 

The cross section for a dark matter particle to scatter 
off of a nucleus with charge Z and atomic number A is 
given by 
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Here \x n is the dark matter-nucleon reduced mass, /„ 
and f p are the relative coupling strengths to protons 
and neutrons, and -Fjv is the nuclear form factor. _Fdm 
is a possible form factor coming from the dark matter 
sector, which is equal to 1 for a standard WIMP, and 
in general is normalized to be 1 at momentum transfer 
q = V2niN-Bj? = 30 MeV. a p is an overall constant that 
would be the dark matter-proton cross section in the case 
where -Fdm = 1- For simplicity we will take f n = in 



this letter, corresponding, for example, to scattering tak- 
ing place through a new C7(l) gauge symmetry which 
mixes with hypercharge. 

The form of the dark matter velocity distribution is 
significantly uncertain. It is generally considered that 
a Boltzmann distribution cx e~ v in the galactic rest 
frame is a reasonable starting point. We will consider 
the cases v = 170km/s, 220km/s, and 270km/s, and take 
the distributions to be cut off at an escape velocity v csc 
whose value we shall take to be 500km/s [T5] . 

We will parameterize the velocity v e of the Earth's 
galactic motion, as well as the nuclear form factors, as in 
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iii. CHANNELING 

Channeling of fast-moving ions through crystal lattices 
is a well-established physical phenomenon with an exten- 
sive literature. For a relatively recent review, see |19j . 
The basic theory, developed by Lindhard in the sixties 
|20j . details under what conditions an ion will become 
guided moving along a particular axis or within a partic- 
ular plane by the cumulative effects of soft collisions with 
the lattice atoms. In this theory, a channeled ion never 
gets close enough to an individual atom to lose a signifi- 
cant amount of energy from elastic scattering. Typically, 
such an ion will instead lose most of its energy by elec- 
tronically exciting the atoms it passes. 

If a WIMP strikes a nucleus in the DAMA detector, 
it is ejected from its position in the Nal lattice (drag- 
ging along most of its electrons [35] ) and becomes a fast- 
moving ion. As pointed out by Drobyshevski [5], this re- 
coiling ion could in principle become channeled, making 
DAMA sensitive to much lower energies than originally 
expected. 

There are two crucial questions here: 1) What is the 
probability for a recoiling ion to become channeled at 
0(1 ~ 10) keV energies in Nal?, and 2) On average, what 
fraction of the recoil energy of such channeled ions is lost 
electronically? In [2T], the DAMA collaboration made 
a first attempt at answering these questions, and using 
Lindhard's theory, estimated the effects to be significant. 
For Na and I ions at 3 keV, they found channeling prob- 
abilities of roughly 30%, and took quenching factors of 
almost 100%. 

However, rigorous experimental and theoretical stud- 
ies of this regime in Nal are lacking, and the full physical 
situation remains far from clear. At DAMA energies, 
calculation of the critical trajectories for channeling falls 
outside the realm of applicability of the Lindhard theory, 
as the details of the collective lattice potential become 
important. Other novel complications appear, as well. 
For example, a recoiling nucleus necessarily originates at 
a lattice site, and thus if it is aimed towards a crystal 
axis or plane, its trajectory will tend to lie in the direc- 
tion of its neighbors. It might then immediately undergo 
scattering away from the intended channel |36j . Even 
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if a recoiling ion manages to find itself channeled after 
a few atomic scatterings with modest energy loss, it is 
not obvious that the subsequent electronic energy losses 
will dominate for such a low-energy channeled ion (see, 
e.g. [12]). Above all else, it is clear that an experimental 
study of the channeling of recoiling nuclei in Nal would 
be extremely valuable. 

In the rest of this paper, we will take an agnostic view- 
point towards these issues, and simply parametrize the 
possibility of a population of recoiling ions depositing an 
anomalously large fraction of its energy electronically. In 
the end, even a somewhat large non-gaussian tail of the 
event-by-event quenching factor can lead to the effects 
predicted by Drobyshevski, and channeling need not even 
be the primary cause of this. Regardless of its origin, we 
will continue to call this hypothetical effect "channeling" 
for continuity. 

Should an effect like this occur at the level claimed 
by DAMA, the implications for their observed annual 
modulation at ~ 3 keV would be quite dramatic. For 
light WIMPs with standard SI interactions, the absolute 
and modulating spectra peak sharply at low energies, 
and DAMA becomes capable of directly sampling part 
of this otherwise unobservable region at Er ~ 3 keV. 
This causes a portion of the DAMA allowed region in 
the (m,(Tjv) plane to edge past constraints from CDMS 
[231 EH ES] and XENON [9]. On the other hand, as 
we will see in the next section, even a modest relaxation 
of either the assumed channeling fraction or its effective 
quenching factor will move all of the allowed region back 
into conflict with these experiments. 

A simple, more robust alternative is a dark matter par- 
ticle with momentum-depedent interactions due to form 
factors or other explicit independence in the scatter- 
ing matrix elements. As we will see below, even modest 
amounts of channeling could make simple form factors in 
the dark sector a viable explanation for all of the present 
data. 



IV. REQUIRED CHANNELING FRACTION 
AND SHAPE 

If the DAMA signal is due to channeled events, then 
the constraints from many other direct detection experi- 
ments can be avoided because most of their energy range 
would correspond to dark matter events with velocity 
above the local escape velocity. It is useful to convert 
the energy range of experiments into a range in momen- 
tum transfers q, with q m j n the lowest momentum transfer 
probed. Then a given experimental constraint is com- 
pletely evaded for dark matter masses below some criti- 
cal value, TOqJ^ = 2(u csc + v e )/q m i n — w^ 1 . In particular, 
it is possible to choose todm so that all of the CDMS-Ge 
and CRESST-W range in q is above the escape velocity, 
while all of the DAMA signal region is below it. The most 
significant constraints then become CDMS-Si and, de- 
pending on its quenching factor, XENON. In addition, it 



/chan~Q-35 




- CDMS-Ge 

- CDMS-Si 

- CRESST-II 

- XENON-lcr 

- XENON 

- ZEPLIN-H 
■■ ZEPLIN-m 

- DAMA unmod 

- COGENT 



m DM (GeV) 

FIG. 1: Constraint plots for channeled events, with an energy- 
independent channeling fraction. The bottom plot shows con- 
straints on models with a momentum-dependent Fum(q) oc q 2 
coupling, and the top assumes a standard (i.e. momentum- 
independent) coupling. DAMA 90% and 99% contours are 
shown. The /chan's above are chosen to be the minimum val- 
ues consistent with all constraints. Obviously, larger values 
of /chan leave more parameter space open. 



is important that values of mrjM below the critical value 
still give a decent fit to the DAMA spectrum. We will see 
that model-dependent assumptions can affect the qual- 
ity of the fit significantly, allowing or disfavoring smaller 
dark matter masses. 

Because the critical mass for a given experiment de- 
pends on g m i n , and event rates tend to rise sharply 
at decreasing q, constraints are very sensitive to even 
modest uncertainties in quenching factors, which relate 
the experimental light yield to the actual recoil en- 
ergy. XENON in particular has measured its quench- 
ing factor to be 15.5 % ± 1.2% (weighted average over 
[26l [27J [28] , but see also [29], which indicates that the ex- 
periments may be underestimating their uncertainties). 
We will therefore present constraints from XENON us- 
ing both (/xenon = 15.5% and the lower la value 14.3% 
("XENON-lcr" for short). 

There are a number of directions in which the pa- 
rameter space opens up if one takes into account un- 
certainties or variations in models. The quality of the fit 
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FIG. 2: Minimum channeling fractions allowed at Er = 
3 keV depending on the channeled quenching factor g c han 
and the channeling fraction energy-dependence, /(Er) oc 
Eft™. Constraints are shown for vo = 170, 220, 270km/s in 
blue dashed, black solid, red dot-dashed, respectively, and 
(/xenon =15.5%(14.3%) in thin (thick) lines. The top panels 
are for a momentum-independent coupling and the bottom 
are for Fdm oc q 2 . Left (right) has a — (g c han = 1.). The 
vo =270km/s constraint in the lower right plot is identical for 
the two values of (/xenon we consider. 



to DAMA depends on the shape of the predicted spec- 
trum, and including an additional dark matter form fac- 
tor Fy)m(<i) ~ Q 2 allows a range of smaller dark matter 
masses. 

It is important to take into account that the masses 
favored by DAMA are also sensitive to the energy- 
dependence of the channeling fraction. Most if not all 
studies of channeled events assume a channeling fraction 
of the form /i(-Er) = f +0 or something very close 
PUDlEn], which over the range 2keV < E R < 6 keV rel- 
evant at DAMA is approximately oc E^ 75 . However, as 
we have discussed, this estimate is based on a relatively 
simple analytic approximation due to Lindhard, and for 
example if the channeling fraction is saturated at the rel- 
evant low recoil energies, then the profile could instead 
be approximately constant |37j . Such a profile would fa- 
vor somewhat smaller masses, which in turn push more 
of the null experiments' energy ranges above the escape 
velocity. Figure [l] shows the 90% constraints on a WIMP 
and the 90%, 99% DAMA fit contours if the channeling 
fraction is energy-independent. A channeling fraction of 
«35% is the minimum value required to have some region 
of parameter space allowed by all constraints. The situ- 
ation is even better with an additional dark matter form 
factor Fdm oc q 2 . In this case, much of the DAMA fa- 
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TABLE I: Minimum channeling fractions at Er — 3 keV al- 
lowed at 90% by all constraints for various values of vo and 
form factors Fdm(?) oc q n . Also shown are the dark matter 
mass at the minimum channeling fraction and, in parenthe- 
ses, the Xdama f° r the best fit to the DAMA spectrum. In 
all cases, (/xenon = 15.5%. 



vored masses are below the critical value ttidm where no 
events are predicted at XENON or CDMS-Ge, and the 
minimum channeling fraction allowed becomes ~12%. In 
all of the allowed parameter space, iodine scattering al- 
ways dominates over sodium scattering at DAMA [35]. 

Our analysis assumptions are as in jl] except for 
XENON, where we vary the quenching factor, CRESST- 
II, where we use only the 2007 run, and CDMS-Si, where 
we have taken an effective exposure of 12 kg days (from 
March 25 - Aug 8) and no events. Also, one of the ten 
potential XENON events was identified by that collab- 
oration as resulting from instrumental error, and we do 
not include it in calculating our constraints. 

Finally, as usual, uncertainty in the halo model has 
a significant impact on the model constraints. We em- 
phasize that it is difficult to know precisely how to in- 
terpret model constraints without better knowledge of 
the allowed range of dark matter halo velocity distribu- 
tions. The most significant effect of different halo models 
tends to be whether the distribution is tighter or broader, 
so one can qualitatively consider the effect of different 
halo models by changing the average rotational veloc- 
ity parameter Vq. In Table [I] we consider the effect of 
a modest change in the halo distribution by considering 
Maxwellian distributions with vo = 170 and 270 km/s 
in addition to the standard vq —220 km/s. We present 
there the minimum channeling fraction at En — 3keV 
allowed by demanding consistency at 90% with all ex- 
periments. Note that, while a higher power of q in the 
form factor decreases the required channeling fraction, 
one may not continue indefinitely to higher powers since 
the best fit to DAMA eventually worsens. Furthermore, 
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when the strongest constraint is XENON, larger values 
of i>o are favored, since the broader velocity distribution 
pushes the DAM A favored region to lower masses. How- 
ever, when CDMS-Si is the strongest constraint, lower 
values of vq are favored, since larger dark matter masses 
enhance the ratio of the reduced mass /i at iodine vs. at 
silicon. 

Figure [2] presents the effect from different possibilities 
for the details of channeling itself, for three halo models 
(v = 170,220,270km/s) and XENON quenching factors 
(feENON = 15.5%, 14.3%). Since it is unlikely that 100% 
of the energy of channeled events gets deposited elec- 
tronically, we consider lower values for this "channeled 
quenching factor" q cna n [39]. For a standard WIMP cou- 
pling, relatively high values (? c han ^ 90% or large chan- 
neling fractions > 30% are necessary. However, for a q 2 
form factor, even g cn an ~70% may be sufficient with as 
little as 15% of events channeled. We also show the effect 
of varying the energy profile of the channeling fraction, 
parameterized as Ji(-Er) oc E~^ . For a standard WIMP 



coupling, the minimum channeling fraction required at 
Er =3keV depends significantly on a, and the profile 
must be fairly flat to allow / c h an <50%. Again, the situa- 
tion is much improved with a form factor. For Fdm °c q 2 , 
the dependence on a is much weaker, and relatively large 

Cd's Still allow /chan ~10%. 
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